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1.  Introduction 

1.1.  In  the  following  we  shall  consider  some  questions  concerning  the  compar¬ 
ison  of  two  samples.  The  test  around  which  our  investigations  will  center  is  the 
Kolmogorov-Smirnov  two-sample  test,  restricted  always  to  the  case  of  equal 
sample  sizes. 

In  the  first  part  we  shall  treat  the  power  function  for  certain  alternatives  and 
make  some  remarks  on  the  efficiency  of  the  test  considered  in  the  case  of  small 
samples.  In  the  second  part  some  remarks  will  be  given  on  distributions  and 
limiting  distributions  occurring  in  connection  with  the  treated  problems.  The 
investigations  given  here  are  closely  connected  with  the  author’s  work  presented 
at  the  Fourth  Berkeley  Symposium. 

1.2.  For  diminishing  the  difference  in  efficiency  between  parametric  and  non- 
parametric  tests,  the  author  has  in  his  papers  [9],  [11]  proposed  the  use  of  a 
pair  of  statistics  instead  of  one  statistic.  In  consequence  of  the  Neyman-Pearson 
lemma,  this  results,  for  given  alternatives,  in  a  better  test  than  the  one  based 
on  either  single  test  statistic.  We  apply  the  maximum  deviation  of  the  two 
empirical  distribution  functions  as  the  first  statistic,  which  ensures  the  asymp¬ 
totic  consistency  of  the  test.  Then  we  can  add  to  this  for  several  types  of  alter¬ 
natives  a  suitable  corresponding  pair,  for  example,  the  first  maximum  index, 
the  number  of  intersections,  the  Galton  statistic,  and  so  on.  In  order  to  examine 
the  increase  in  the  efficiency  of  the  two-sample  Smirnov  test,  we  shall  treat  the 
situation  in  the  case  of  a  special  alternative,  for  which  the  computation  is 
relatively  easy. 

In  our  treatment  we  make  use  of  the  power  functions  of  the  original  test  and 
of  the  two-statistic  test  as  well.  The  power  function  can  be  constructed  easily 
in  case  of  a  (continuous)  alternative  containing  piecewise  linear  parts.  With 
such  alternatives  we  can  approximate  any  given  alternative.  Following  Z.  W. 
Birnbaum  [1],  these  kinds  of  alternatives  (for  instance,  the  maximum  and  min¬ 
imum  alternatives)  were  treated  by  many  authors  in  the  one-sample  case.  As 
we  shall  see  in  section  1,  this  power  function  can  be  easily  obtained  in  the  two- 
sample  case  for  all  tests  for  which  the  distribution  of  the  test  statistic  under 
null  hypothesis  is  known;  the  idea  used  is  the  extension  of  the  method  used  by 
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Z.  W.  Birnbaum  in  his  mentioned  paper.  Nevertheless,  this  is  simple  enough 
for  further  considerations  only  when  the  number  of  linear  sections  of  the  dis¬ 
tribution  function  is  small.  Further,  we  need  the  power  function  for  the  more 
complicated  case  of  two-statistics  tests,  so  we  shall  work  in  this  first  occasion 
with  a  very  simple  alternative  knowing  that  our  aim  is  to  discover  what  we  can 
expect  at  all  from  the  use  of  pairs  of  statistics.  We  shall  consider  the  power 
for  moderate  or  small  sample  sizes  which  occur  very  often  in  the  applications. 
Our  formulae  give  the  possibility  of  carrying  out  the  program  of  D.  Chapman  [2] 
which  was  done  for  the  one-sample  case. 


2.  The  power  of  two-sample  tests  for  piecewise  linear  alternatives 

2.1.  Notations.  Let  ^  and  tj  be  random  variables  with  continuous  distribution 
functions  F{x)  and  G{x)  resp.  for  which  the  null  hypothesis  Ho:  G{x)  and  F{x) 
are  equally  and  uniformly  distributed  in  interval  (0,  1).  We  shall  treat  piecewise 
linear  alternatives  depending  on  vectors: 


(2.1) 

Z  =  (zo,  2l,  •  •  •  ,  2r), 

(0  =  Zo  <  2l  <  22  <  •  • 

•  <  Zr  =  1), 

(2.2) 

g  =  (gi,  gi,  , gr), 

(gi  >  0,  i 

T 

:  =  1,  2,  •  •  •  ,  r)  with  E  gi(zi 

i=l 

^i—j)  ~  f  • 

C F(x)  as  in  Ho, 

(2.3) 

|G'(a:)  =  g(x)  =  - 

fg,  if  Zi-i  <  X  <  Zi,i  =  1,  2,  • 
1.0  otherwise. 

•  •  ,  r. 

In  this  case  g{x)  dx  =  gix)  dx  =  1. 

Let  further  ^i,  I2,  •  •  *  ,  and  rji,  ,  Vn  be  independent  observations  on  ^ 

and  77  resp.  We  denote  the  elements  of  the  ordered  samples  by  and  77*  resp., 
and  the  union  of  the  two  samples  in  order  of  magnitude  by 

(2.4) 

T*l  < 

Tt<  <  Tin. 

Let  us  define  further  the  random  variables  for  f  =  1,2,  ■  •  •  ,  2n, 

(2.6) 

di  =  j 

'+1,  if 

—  1,  if  =  77/ 

for  some  j  and  With  the  usual  notation  so  =  0,  Si  =  +  •  •  •  + 

(i  =  1,  2,  •  •  •  ,  2n),  S2n  =  0.  The  points  (i,  Si)  in  the  plane  give  the  path  of  a 
random  walk  starting  at  the  origin  and  returning  after  2n  steps  to  the  point 
(2n,  0). 

In  the  present  paper  we  shall  consider  the  following  statistics : 

(2.6)  =  max  {Fn{x)  —  G„{x))  =  -  max  s,-, 

(x)  n  (t) 

(2.7)  Dn,n  =  max  |F„(x)  —  (j„(a:)|  =  -  max  \si\, 

(x)  ^  (i) 
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(2.8)  =  Si  =  nD^nj* 

(2.9)  A„.n  =  the  number  of  intersection  points  in  the  above  mentioned  random 

path,  that  is,  the  number  of  i’sfor  which  Si  =  0  and  Si-iSi+i  =  —  1 
occurs,  adding  the  point  (2n,  0). 

As  to  the  numerical  determination  of  the  values  of  these  statistics  we  mention 
the  following.  As  it  was  pointed  out  in  [3],  nD^n  agrees  with  that  index,  for 
which  in  the  translation  scheme 


(2.10) 


•  •  •  >  k*,  kt+i,  *  •  •  ,  k*+s,  '  "  ,  k 


* 

n 


Vi, 


Vs, 


Vn-» 


Vn 


it  first  occurs  that  each  ^*+i  exceeds  the  corresponding  tj*.  Translation  in  the 
opposite  direction  leads  to  D^,„  and  in  this  way  to  D„,n.  Further  =  k  +  2s 
where  s  =  min  {i:  <  17*+ 1},  which  can  be  seen  easily. 

Having  the  two  ordered  samples,  let  us  define  the  random  variables 


(2.11) 


r+i,  if  ^1  >  17*, 

\-i,  if  ei<vi 


Then  A„,„  —  1  equals  the  number  of  changes  of  sign  in  the  sequence  ei,  62, 
•  •  •  ,  6„  (see  [12]). 

2.2.  The  power  function  of  the  two-sample  Smirnov  test  under  the  hypothesis 
Hq  against  H\''.  As  is  known  the  a-size  critical  region  for  one-sided  alternatives 
F{x)  =  G{x) — this  case  will  be  treated  in  more  detail  later  in  this  article — is 
determined  by  the  relation 


(2.12) 


a 


with  k  =  k{a,n). 

Let  us  denote  by  the  vectors  v  =  {vi,  vz,  •  ”  ,  Vr)  and  p  =  (m,  mz,  •  •  *  ,  Mr) 
the  events  that  out  of  the  n  exactly  vj  and  (independently)  out  of  the  n  rjiS 
exactly  m  are  contained  in  the  intervals  (zy-i,  zj)  j  =  1,  2,  ,  r.  The  prob¬ 

abilities  of  these  events  are  clearly 


P{v\Ho)  =  P{v\Hi)  =  n!  n 


(Zj  -  Zy-i)"’- 


(2.13) 


j  =  i  vfi 


=  nth  — - 

3=1  Md 


(2.14)  PMH'P)  =  n!  ri  ~ 

3  =  1  P’3  • 

further,  for  f  =  0,  1, 


(2.15)  P{v,  p\Hi)  =  P{pmP{p\Hi). 


Denoting  now  for  j  =  1,  2,  •  •  •  ,  r  by  wy  and  wy  the  partial  sums  j'i  +  i'2  + 
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•  •  •  Pj  and  Ml  +  Ms  +  •  •  *  +  Mi  resp.,  (no  =  mo  =  0),  then  the  following  holds. 

Theorem  2.1.  For  the  a-size  one-sided  two-sample  Smirnov  test  the  power 
^unction  under  Hi^  against  Ho  is  given  by 

(2.16)  WniHV,  a)  =  1  -  (n!)2  E*  E*  H 

(f)  (ji)  i = 1 

/  +  Mi 

\My  +  fc  —  ny_i  +  my-i 

where  Ew  E(m)  denote  r-summaiion  for  all  possible  v  and  /i  vectors  with  0  <  Vj, 
Mi  <  n  and  =  Ei=i  Mi  =  ^  the  further  restriction  nj  —  m,  <  k,  j  = 

1,  2,  •  •  •  ,  r. 

Before  outlining  the  proof  we  make  some  remarks. 

(a)  The  expression  of  the  power  function  seems  to  be  suitable  for  asymptotic 
considerations  or  computational  work  in  general  for  small  values  of  r.  For  r 
large  and  for  small  intervals  (2,_i,  2,)  we  are  interested  in  the  mutual  order  of 
sample  elements  inside  the  intervals  only  for  very  peculiar  alternatives.  In  these 
cases  the  problem  can  be  reduced  into  the  consideration  of  the  null  hypothesis 

(2.17)  Hh:  PU  e  (zi-i,  Zi))  =  Pin  e  (2,_i,  2,))  for  i  =  1,  2,  •  •  •  ,  r. 

This  means  that  our  problem  is  the  comparison  of  two  multinomial  distributions, 
which  is  treated  recently  by  Hoeffding  [6]. 

(b)  If  n  remains  finite  and  r  tends  to  infinity,  then  we  come  essentially  to 
the  evaluation  of  the  probabilities  of  each  different  array  of  the  sample  elements, 
arrays  which  belong  to  the  critical  region.  This  kind  of  expression  was  given 
by  Hoeffding  [5]  and  considered  by  Lehmann  [7]. 

(c)  In  the  two-sided  case  the  combinatorial  quantity  just  after  the  product 
sign  is  to  be  replaced  by  the  expression  corresponding  to  the  two-barrier  case 

“  {vi  +  -  mti  +  (27  +  m 

where  for  the  summation  the  restriction  \nj  —  rrijl  <  k,  j  =  1,  2,  •  •  •  ,  r  is  to  be 
made. 

The  proof  is  the  consequence  of  some  simple  arguments  which  are  used  often 
for  similar  purposes  and  which  are  the  following. 

Lemma  2.1.  Let  be  a  continuous  random  variable  in  the  interval  (0,  1)  whose 
density  function  is  constant  in  a  subinterval  (a,  b).  Then  is  uniformly  distributed 
in  (a,  b)  under  the  condition  {a  <  ^'  <  b}. 

A  consequence  of  this  is  the  following  lemma. 

Lemma  2.2.  Let  and  n'  be  continuous  random  variables  in  the  interval  (0,  1), 
the  density  functions  of  which  are  {not  necessarily  equal)  constants  in  a  subinterval 
(a,  6).  Let  ^1,  ^2,  •  •  •  ,  and  tjI,  172,  •  •  •  ,  t/m  be  independent  observations  on  and 


)(zj  - 
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rj'  resp.  falling  in  the  interval  (a,  6).  Under  this  condition,  all  possible  arrays  of 

the  mentioned  v  -{■  n  sample  elements  have  the  common  probability 

A  consequence  of  this  lemma  is  that  in  each  interval  the  conditional  prob¬ 
ability  that  Si  <  k  given  {/»,  v} ,  can  be  calculated  under  Hq. 

Turning  now  to  the  determination  of  the  probability 

(2.19)  1  -  a)  =  P  ^  IHf’V 

we  calculate  this  conditionally  given  {/t,  p}  and  multiply  this  by  P{p, 
then  we  have  to  sum  over  all  v  and  /i.  But  under  the  mentioned  condition  the 
random  walk  falls  into  parts  with  division  points  {ij  =  w,-  -f-  my,  sy,  =  ny  —  my), 
{j  =  1,  2,  •  •  •  ,  r)  and  according  to  the  Markov  property  of  the  random  walk, 
the  relation 

(2.20)  P  (ci.  <  ^  1  M,  »,  Hi'') 


=  n  P(s*  <  k,  for  ny_i  -}-  my_i  <i  <  nj  my|/i,  v,  Hi^) 
i=i 

holds.  Taking  into  account  the  elementary  formula 
(2.21)  P{si  <  k,  for  wy_i  4-  my_i  <  i  <  nj  -f-  myl/t,  v, 


(  Vj  -I-  Mi  \ 

—  1  _  \Mi  +  fe  ~  nj-i  -f-  my-i/ 


and  the  relations  (2.13),  (2.14),  and  (2.15),  we  come  to  theorem  2.1. 

2.3.  A  special  case.  For  our  comparative  considerations  we  shall  treat  the 

uniformly  distributed  in  (1, 1), 

0  if  0  <  X  <  s, 

,  I Y— —  if  z  <  x  <  1. 

L  LI  —  z  ~ 

For  this  kind  of  alternative  the  first  index  of  the  maximum  will  occur  with 
higher  probability  for  smaller  values  than  in  the  case  of  the  null  hypothesis. 
This  alternative  being  simple  enough,  I  have  chosen  it  for  a  first  comparison 
of  tests  based  on  1  and  on  2  statistics  respectively. 

Making  use  of  the  notation  vi  =  v,  then  V2  =  n  —  v,  and  further,  knowing 
that  P(mi  =  0,  M2  =  n\H^^)  =  1,  we  obtain  for  the  power  function 


following  very  simple  alternative 
(Fix) 


(2.22)  ] 


nl(rr\  —  _ 


2n  —  V 
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This  simple  case  shows  that  in  determining  the  power,  the  normal  approx¬ 
imation  of  the  binomial  terms  is  not  suitable,  since  v  <k{a,n)  ^  when 

n  — >  00 .  A  similar  remark  was  made  by  J.  Rosenblatt  [8]. 

After  some  modification  we  may  obtain  the  following  form: 


(2.24)  WM"?,  a)  =  1  -  Bi-=(«  +  1,  *:), 

where  B  is  the  beta-function 


(2.25) 


Bziv,  q) 


—  0®“^  dt 


This  is  suitable  for  inunediate  computation  of  the  power  against  this  simple 
kind  of  alternative.  Now  we  give  some  numerical  values  for  small  samples  using 
(n,  k)  pairs  for  which  a  is  near  10%  and  z  =  max(x)  {F{x)  —  G{x))  =  0.1; 
0.2;  0.3.  The  following  table  gives  the  errors  of  the  second  kind  for  these  values 
of  2. 


TABLE  I 


n 

k 

a 

z  =  0.1 

z  =  0.2 

II 

o 

20 

7 

0.0873 

0.8177 

0.6163 

0.3363 

30 

8 

0.1197 

0.7230 

0.4039 

0.1164 

40 

9 

0.1331 

0.6587 

0.2640 

0.0377 

50 

11 

0.7190 

0.7190 

0.2700 

0.0200 

2.4.  The  power  function  of  the  test  based  on  a  pair  of  statistics.  Let  us  consider 
the  two  statistics  Dn,n  and  Rf,n  for  the  decision  between  Ho  and  defined  in 
2.3.  We  shall  introduce  the  random  variable  =  Rf^,n  —  Dt,n.  As  nD^»  and 
nRt,n  are  of  the  same  parity,  Stn  is  always  even.  Let  us  use  the  following  nota¬ 
tions  : 

(2.26)  P  (bJ.  =  S.+.  =  ^  I  Ho)  =  PiS 

and 

(2.27)  P  (bJ.  =  St.n  =  ^  1  Bi)  =  Qi5. 

Denoting  the  best  critical  region  of  the  a-size  test  (restricting  ourselves  to 
the  {k,  s)  plane)  by  3Ca,  this  is  defined  with  the  aid  of  a  suitable  constant  Ca, 
and  can  be  written  in  the  form 


(2.28) 


3Za  = 


>01 

p,„>  >  c.j- 


and 

(2.29) 
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E  ra  = 

(*!,«)  eac- 


The  power  function  is 

(2.30)  Wn{H?\  a)  =  E 

(A:,«)e3C» 

We  need  the  probabilities  P^l  and 
As  it  was  proved  in  [10]  the  following  relations  are  valid: 


(2.31)  P67  = 


(2.32)  P’ts  = 


1  1 

f'2s\ 

Ks) 

/2n 
\  n 

2(2s  -  l)(n  -  s  +  1) 

/2n\ 
\n } 

1 

s  =  1,  2,  •  •  •  ,  n, 


k{k  +  1) 


(‘r)( 


—k  +  2w  —  2s 
n  —  s 


(k  +  2s)(n  —  s  +  1) 


), 


s  =  0,  1,  2,  •  •  •  ,  n  —  A:. 

We  turn  now  to  the  determination  of  the  probabilities  The  number  of  ^/s 
in  the  interval  (0,  z)  may  be  v  =  0,  1,  2,  •  •  •  ,  A:.  As  was  mentioned,  /ii  =  0  with 
probability  1.  For  A:  =  0,  j'  must  be  0  and  P{v)  =  P{v  =  0,  n  —  v  =  n)  = 
(1  —  2)".  In  this  case  it  can  be  seen  inmiediately  that 


(2.33) 


«  =  p^ri(i  -  zy, 


s  =  1,  2,  •  •  •  ,  n. 


The  case  s  =  0  (that  is,  Rn,n  =  Dn.n)  can  happen  only  when  v  =  A:.  In  this 
case  we  have  the  second  part  of  the  path  starting  with  probability  1  from  the 
point  {k,  Sk  =  k)  and  ending  at  (2n,  0).  The  probability  that  this  path  will  never 

reach  the  height  s,-  =  A:  +  1  multiplied  by  P{v)  =  P{v  =  k)  =  2*(1  —  2)”“* 

gives  the  required  probability 


(2.34) 


-  2)"-* 


A;  =  1,  2,  •  •  •  ,  w. 


At  the  end,  for  s  >  0,  A:  >  0,  we  can  construct  the  power  function  as  given 
in  2.2  and  evaluate  the  joint  probabilities  for  the  maximum  and  first  maximum 
index  in  the  case  of  the  several  p’s. 
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The  resulting  formula  is 

.  (' 


(2.35)  = 


k  + 


■fc  +  2w  —  2s 
n  —  s 


)  k  —  V 


w  —  s  +  1 


© 


=0  n 


Using  the  above  formulas  for  numerical  calculation  we  can  compare  the 
second  kind  of  errors  for  the  use  of  Bt.n  alone  and  of  the  pair  ipt.n,  Rn,v)  respec¬ 
tively.  The  results  are  tabulated  in  table  II.  (The  computation  was  carried  out 
on  a  GIER  electronic  computer.  I  am  indebted  to  A.  B4k6ssy  and  G.  Tusnddy 
for  their  kind  help  in  accomplishing  these  calculations.) 


TABLE  II 


n 

Error  of  first  kind 

Dn,n  Rn.n) 

Error  of  second  kind 

Dn,n  {Dn,nf  Rn.n) 

if  A  =  Z  =  0,  1 

[  in  the  case  of  using 

Dn.n  {Dn.tit  Rn[n) 

if  A  =  z  =  0,  2 

10 

0.0739 

0.0839 

0.8581 

0.7532 

0.6198 

0.2005 

0.2038 

0.6815 

0.5215 

0.3683 

30 

0.0675 

0.0675 

0.0893 

0.8263 

0.6149 

0.5636 

0.5573 

0.2568 

0.1197 

0.1197 

0.7230 

0.5036 

0.4039 

0.1703 

50 

0.0562 

0.0562 

•0.8017 

0.5105 

0.3812 

0.1179 

0.0893 

0.0893 

0.7190 

0.4299 

0.2700 

0.0780 

0.1362 

0.1362 

0.6133 

0.3488 

0.1733 

0.0473 

3.  The  mflyimum  and  the  number  of  intersections 

3.1.  The  nonconsistency  of  the  number  of  intersections.  In  our  paper  with  E. 
Csdki  [3]  we  considered  the  statistic  A„,„  (see  2.1),  that  is,  the  number  of  inter¬ 
sections.  As  a  test  statistic  this  has  the  following  advantages.  As  we  mentioned 
in  1.1,  its  value  can  be  determined  very  easily.  The  distribution  of  A„,„  under 
Hq  is  very  simple  too : 


In  addition  to  these,  the  test  based  on  A„,„  has  the  same  properties  for  one-sided 
and  two-sided  alternatives. 

We  conjectured  that  this  statistic  is  consistent  against  all  continuous  alter¬ 
natives.  The  grounds  for  this  conjecture  were  the  following :  if  max(x)  {F{x)  —  G{x)) 
=  A  >  0,  then  with  probability  greater  than  zero  the  statistic  will  take 
values  of  order  of  magnitude  n.  Now  the  following  theorem  is  valid. 
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Theorem  3.1  (E.  Csdki  and  I.  Vincze).  Let  di,  ^2,  •  •  •  he  independent  ran¬ 
dom  variables  with  P(di  =  1)  =  P(^,-  =  —  1)  =  1/2.  Let  So  —  0,  Si  =  + 

^2  +  '  •  •  di,  i  =  1,2,  •  •  •  ,  and  let  us  define  X„  as  the  number  of  i's  for  which 
(Si  =  0,  (Si_i(Si+i  =  —  1,  1  <  i  <  n.  Then  the  following  relation  holds  with 
0  <  c  <  1, 

(3.2)  lim  P(X.  =  <|S„~cn)  =  (  =  0,1,2,  ■  ■  ■. 

Further  if  }f({n)  — > «  and  <  1,  then 

(3.3)  lim  P  ^X„  <  y\Sn  ~  =  1  -  y  >  0. 

Now  from  this  argument  we  would  think  that  for  n  large  enough  the  case 
Fix)  =  Gix)  will  lead  to  the  greatest  number  of  intersections  (in  probability), 
and  thus  the  critical  region  will  be  the  small  values  of  An.n- 
But  E.  M.  Sarhan  (unpublished)  has  given  an  example  that  A„,n  is  not  con¬ 
sistent  against  the  following  alternative  in  (0, 1)  for  2  >  0, 

(3.4)  Fix)  =  X,  0  <  a;  <  1 


(3.5) 


fx, 

Gix)  =  j  -f  z 

llTi’ 


if  0  <  X  <  2, 
if  2  ;<  X  <  1. 


On  the  other  hand,  he  showed  that  the  test  based  on  A„,n,  using  for  the  decision 
between  Ho  and  defined  in  2.3,  is  more  efl&cient  than  the  one-sided 
Kolmogorov-Smirnov  test.  This  way  An.n  as  a  test  statistic — by  itself  or  in 
addition  to  the  Smirnov  statistic — seems  not  to  be  without  interest. 

3.2.  Joint  distribution  of  the  maximum  deviation  and  the  number  of  intersections. 
In  our  paper  [4]  with  E.  Csaki  the  generating  function  of  Pn.n  and  An.n  is 
determined  under  Ho,  which  is  the  following 


(3.6) 


where 

(3.7) 


1  -  Vl  -  42 
^  ~  1  +  Vl  -i-  42* 


We  can  obtain  without  any  difficulty  the  probabilities  by  series-expansion; 
these  are  the  following: 

j  — 


(3.8) 


P  {d„,„  <  A.,„  =  (\Ho)  =  7^  f  t 

\  n  /  / 2w\  j—o 


i-iy 


ro 


X 


'  2n  ^ 

\ji  +  iifi  —  1)  +  jik  -j-  1)  + 


^•(fc-l)-^j(fc-|-l)+^ 

n 
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This  formula  is  not  suitable  for  the  determination  of  the  limiting  distribution 
as  n  00 ,  because  in  the  sum  each  term  tends  to  infinity.  The  way  of  solving 
the  problem  was  the  evaluation  out  of  the  integral  form 


where  the  integration  path  is  a  small  circle  around  the  origin.  This  was  kindly 

done  by  N.  G.  de  Bruijn,  which  we  give  in  the  following  theorem. 

Theorem  3.2  (N.  G.  de  Bruijn).  If  x  >  Q,  y  >  Q,  then 

(3.10)  lim  P  Dn,n  <  y,x  <  A„.n  <  X 

=  i-1  "’'P  ^ 

hold. 

To  the  proofs  of  theorems  3.1  and  3.2  and  a  detailed  treatment  of  the  ques¬ 
tions  of  the  joint  distribution  law,  we  should  like  to  return  later. 
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